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Abstract. A geometric analysis of the sdg interacting boson model is performed. 
A coherent-state is used in terms of three types of deformation: axial quadrupole 
(fa), axial hexadecapole and triaxial (72)- The phase-transitional structure 

is established for a schematic sdg hamiltonian which is intermediate between four 
dynamical symmetries of U(15), namely the spherical U(5) ® U(9), the (prolate and 
oblate) deformed SU±(3) and the 72-soft SO(15) limits. For realistic choices of the 
hamiltonian parameters the resulting phase diagram has properties close to what is 
obtained in the sd version of the model and, in particular, no transition towards a 
stable triaxial shape is found. 



1. Introduction 

Phase transitions are well-known phenomena. They are said to occur if an order 
parameter of the system experiences, as function of a control parameter, a rapid change 
at a particular point. A phase transition may be of first or second order, depending 
on the discontinuity of the first- or second-order derivative of the order parameter as 
a function of the control parameter (Ehrenfest classification). Phase transitions have 
been observed in countless classical and quantal systems, usually composed of very large 
numbers of particles and hence allowing a statistical treatment. The atomic nucleus, 
however, consists of less than 300 particles (nucleons) and finite-system effects may 
therefore play a role. 

The concept of phase transitions is used in nuclear physics either to describe the 
behaviour of hot nuclei resulting from heavy-ion collisions or in the context of so-called 
'quantum' phase transitions, also called shape phase transitions. The latter is the subject 
this paper and relies on a method proposed by Gilmore in 1979 pQ. Its central idea is that 
any quantum-mechanical hamiltonian that can be written in terms of the generators of a 
compact Lie algebra, has an unambiguously defined classical analogue with a geometric 
interpretation in terms of shape variables of which the phase-transitional behaviour can 
be studied. A few years before Gilmore's proposal, in the middle of the 1970s, Arima 
and Iachello [2j [3j HI [5] had proposed a novel way to describe the nucleus with a set of 
s and d bosons, called the interacting boson model (IBM) [6], which precisely has the 
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property that it can be cast in terms of generators of the unitary Lie algebra U(6) to 
which Gilmore's technique can be applied. This approach was advocated by Dieperink et 
al. [7j and lead to a geometric interpretation of the IBM [HI [9] and to a connection with 
the Bohr-Mottelson description of the nucleus [10]. The problem of phase transitions 
in the IBM was subsequently addressed by Feng et al. [llj but was solved in its full 
generality only many years later by Lopez- Moreno and Castanos |12j . 

Although the early work in the 1980s determined the most essential phase- 
transitional properties of the IBM, many additional related features of the model have 
been studied over the years [131 HH US] • This area of research received renewed impetus 
following Iachello's suggestion of critical-point symmetries E(5) and X(5) [161 Eji 
nuclear examples of which were proposed by Casten and Zamfir [TBI [19] . The central idea 
of Iachello's method is to provide a description of nuclei at the phase-transitional point 
through the (approximate) analytic solution of a geometric Bohr hamiltonian. This 
rekindled interest in the phase-transitional behaviour of algebraic models from various 
angles: the correspondence between the E(5) solution of the Bohr hamiltonian and the 
U(5)-SO(6) transition in the IBM was studied in detail [201 1211 [221 [231 121] , the existence 
of an additional prolate-oblate transition was recognized [251 [26] . a connection with 
Landau theory of phase transitions was established [23 [2S [22] , the influence of angular 
momentum was investigated [301 [311 [32], the connection with shape coexistence was 
revisited [33], the algebraic framework of quasi-dynamical symmetries was established 
to deal with transitions between two phases [3U [351 ESI [371 EH], phase transitions of 
excited states were considered [39]. In addition, several extensions of the IBM were 
re-examined from this angle such as the neutron-proton version of the model [10] or its 
configuration-mixed version [4Tj . 

In this paper we report on a phase-transitional study of another extension of the 
IBM, namely the sdg-IBM, which includes a g boson to account for hexadecapole 
deformations of the nucleus. We give a brief review of the sdg-IBM in section [2] and 
discuss its classical limit in section El Since a general sdg hamiltonian has too many 
interaction parameters, a simpler parametrization is proposed in section HI The method 
for establishing the phase diagram is explained in section [5] for the specific case of the 
sdg-IBM and then applied to obtain partial and complete diagrams in sections [6] and [7j 
respectively. The conclusions of this work are summarized in section [HI 

2. The general hamiltonian of sdg-IBM 

We refer the reader to the paper of Devi and Kota 02] for an excellent review of studies 
with the sdg-IBM. We limit ourselves here to a concise description of the hamiltonian 
of the model which is needed in this work. The building blocks of the sdg-IBM are s, 
d and g bosons with angular momenta £ = 0, 2 and 4. A nucleus is characterized by a 
constant total number of bosons N which equals half the number of valence nucleons 
(particles or holes, whichever is smaller). No distinction is made here between neutron 
and proton bosons. 
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Since the hamiltonian of the sdg-TBM conserves the total number of bosons, it can 
be written in terms of the 225 operators b\ m bi' m i where b\ m (b em ) creates (annihilates) a 
boson with angular momentum I and z projection m. This set of 225 operators generates 
the Lie algebra U(15) of unitary transformations in fifteen dimensions. A hamiltonian 
that conserves the total number of bosons is of the generic form 

H = E + H {1) + H {2) + ■■■, (1) 

where the index refers to the order of the interaction in the generators of U(15). The 
first term Eq is a constant which represents the binding energy of the core. The second 
term is the one-body part 

HW = Cs [ s t x 5 ](o) + Qv /5[rft x J](o) + €g yfgy x g\i°) 

= e s n s + e s n d + e d n g , (2) 

where x refers to coupling in angular momentum (shown as an upperscript in round 
brackets), b em = [—)^ m b t _ m and the coefficients e s , e d and e g are the energies of the 
s, d and g bosons, respectively. The third term in the hamiltonian (TjQ) represents the 
two-body interaction 

H (2) = E >< b j J L) x [b t , x \]^\ (3) 

where the coefficients v are related to the interaction matrix elements between 
normalized two-boson states, 



or _i_ i 

'\i 2 -LM\E^\e x e 2 -LM). (4) 



Since the bosons are necessarily symmetrically coupled, allowed two-boson states are 
s 2 (L = 0), sd (L = 2), sg (L = 4), d 2 (L = 0,2,4), dg (L = 2,3,4,5,6) and g 2 
(L = 0, 2, 4, 6, 8). Since for n states with a given angular momentum one has n(n + l)/2 
interactions, 32 independent two-body interactions v are found: six for L = 0, ten for 
L = 2, one for L = 3, ten for L = 4, one for L = 5, three for L = 6 and one for 
L = 8. Together with the constant E and the three single-boson energies e s , e d and 
e g , a total of 36 parameters is thus needed to specify completely the hamiltonian of the 
sdg-IBM which includes up to two-body interactions. This number is far too great for 
practical applications. In the following sections we explain how possible simplifications 
draw inspiration from the classical limit and from the s<i-IBM. 



3. The classical limit 



The classical limit of an IBM hamiltonian is defined as its expectation value in a coherent 
state p] and yields a function of the relevant deformation parameters which can be 
interpreted as a potential surface depending on these parameters. The method was first 
proposed for the s<i-IBM [HJ [7J |9]. The extension to the sdg-IBM was carried out by 
Devi and Kota [43] who established the classical limit of the various limits of U(15). 
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The coherent state for the sdg-IBM is given by 



N 



N; a 2M , a 4M ) oc s 1 " + ^ a 2 ^ + ^ a ^gl |o), (5) 



where |o) is the boson vacuum and the are the spherical tensors of quadrupole 
(A = 2) and hexadecapole (A = 4) deformation. They have the interpretation of 
shape variables and their associated conjugate momenta. If one limits oneself to static 
problems, the a>\^ can be taken as real; they specify a shape and are analogous to the 
shape variables of the droplet model of the nucleus [10]. In particular, they satisfy 
the same constraints as the latter which follow from reflection symmetry, that is, (i) 
= and (ii) = for \i odd. 

The most general parametrization of a surface with quadrupole and hexadecapole 
deformation has been discussed by Rohozihski and Sobiczewski jH] and involves five 
shape parameters: two of quadrupole (f3 2 and 72) and three of hexadecapole character 
(/?4, 74 and £4). A considerable simplification of the problem is obtained by relating the 
two hexadecapole variables (74 and 64) that parametrize deviations from axial symmetry 
in terms of the corresponding quadrupole variable 72, according to the Cayley-Hamilton 
theorem [35]. In terms of the latter variables, the coherent state ([5]) is rewritten as 



2 -T M4 ; 

t 



x (s + cos 72^0 + V \ sin 72(^2 + d\ 2) 

+ \(5 A [(5 cos 2 72 + l)g\ + M sin 2 l2 {gl 2 + g\ 2 ) 



+ v /f sin 2 7 2 (£ 4 + ^4)Jj"». (6) 

The range of allowed values for the quadrupole shape variables is < (3 2 < +°o 
and < 72 < tt/3. On the other hand, the hexadecapole shape variable /3 4 is 
unrestricted [15]. The expectation value of the hamiltonian ([1]) in this state can 
be determined by elementary methods |46j and yields a function of /3 2 , Ai and 72 
which is identified with a potential V(02, 04, 72)- In this way the classical limit of 
the hamiltonian (CQ) has the following generic form: 

E<L(ft) l (A)'cos(3m 72 ), (7) 



n>l 
X 



fcZm 

where n refers to the order of the interaction. If up to two-body terms are included in 

An) 
L kl;m 



the hamiltonian, the non-zero coefficients 



(i) a) (i) 

a 00;0 — £ si a 20;0 — e di <3- 02; o — c g , 

a 00;0 _ 2 U ssssi a 20;0 _ y 5 u ssdd ' ^sdsd' a 02;0 _ 3 ^ssgg + ^sgsff' 

„( 2 ) _ O /T„.2 „(2) _ f2„,2 

"30:1 



z y 7 U sdddi a 21;0 ~ z y 7 y s(i<2c/ + y /^ u sgdd 



(2) _ _10 /T 2 _ 4 /T„4 „(2) _ 112 ,4 

l 12;l - 3,\l77 u sdgg ^\l77 U sgdg-> «03;0 ~ 9^/1001 
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(2) _ 50 4 (2) _ i n , 12 , _9_ 4 

"03;2 — gyTooT 40 ;° ~~ 10 U dddd ' -j u dddd * ^ u ddddi 

(2) _ _2s/2 2 12_ 4 

a 31;l _ 7 u dddg 7 ^/fi U dddg^ 

a 22;0 - 3v ^ V ddgg + 9 \ U V ddgg + 3VnE V ddgg 

_L 2 2 J_ 3 _85_ 4 , 47 6 

~ 7 rf9fi9 ~ 18 ^9 i 462 ' 99 u dgdg 5 

„( 2 ) = J. 4 /X,,a + 1° /IT?, 4 _ J_ ? ,3 , _5_,,4 _ 2.,.6 
"2212 63 V H ^99 ~ 21 V 143 u ddgg 18 u dgdg ~ m dgdg 99 v dgdg i 

J2) _ 20 „,2 36 /~5~„,4 10>/2„,6 



t13 ; 1 2lVn V dggg 77 \ 13 V dggg 33 ^999' 

(2) = J_ 1 _625_ ? ,2 I 167 4 , _5_ 7 ,6 1 2485 8 

"04;0 18 9999 ~ 6237 9999 ~ 2457 9999 ~ 81 9999 ~ 11583 9999' 

(2) 25 2 I _50_ ? ,4 I _35_ 7 ,6 , 280 8 /o\ 

u 04;2 891 L 9999 ^ 3861 U 9999 ^ 891 ^9999 11583 U 9999' 

in terms of the single-boson energies q and the matrix elements between normalized 
two-boson states, 



2 



LM\H {2) \t'^ 2 -LM). (9) 



This analysis shows that the geometry of the sdg hamiltonian depends on 19 
parameters (including Eq), considerably smaller than the 36 parameters of the quantum- 
mechanical version of the same hamiltonian. The structure of the phase diagram 
associated with the general sdg hamiltonian ([7]) might be of interest but is well beyond 
the scope of the present analysis. In the next section we show how a simplified sdg 
hamiltonian can be constructed which presumably captures the essential features of the 
phase-transitional behaviour in the sdg-lBM. 



4. A generalized Casten triangle 

A simplified hamiltonian of the sd-IBM is of the form 

H sd = e d h d + kQ ■ Q + k'L ■ L, (10) 

where Q is the quadrupole operator with components = [^xs+s^x <^]L 2 ^+x[<^ x d]^ 
and L is the angular momentum operator, = vlo [S x d]^\ The Q 2 and L? terms 
in fflOl) constitute the hamiltonian of the so-called consistent-Q formalism (CQF) [47J; 
for x = ±V7/2 it gives rise to the deformed or SU(3) limit and for x = to the 
7-unstable or SO(6) limit. In an extended consistent-Q formalism (ECQF) [48J the 
term fid is added with which the third, vibrational or U(5) limit of the sd-IBM can be 
obtained. The ECQF hamiltonian thus allows one to reach all three limits of the model 
with four parameters. 

The eigenfunctions of the ECQF hamiltonian are, in fact, independent of k' and, 
furthermore, of the overall scale of the spectrum, reducing the number of essential 
parameters to two, the ratio K/ea and \- The physically relevant portion of this 
parameter space can be represented on a so-called Casten triangle [49] . Its vertices 
correspond to the three limits of the sd-IBM, U(5), SU(3) and SO(6) for spherical, 
deformed and 7-unstable nuclei, respectively. The three legs of the triangle describe 
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transitional cases, intermediate between two of the three limits, while an arbitrary point 
on the triangle corresponds to an admixture of the three limits. 

In this paper we propose a similar hamiltonian in the sdg-IBM. The different 
dynamical symmetries of the sdg-IBM are known since long [50l |5T] . In particular, three 
limits occur in sdg-IBM, namely U(5)®U(9), SU(3) and SO(15), that are the analogues 
of the limits of the sd-lBM. An immediate generalization of the sd hamiltonian fflOl) is 
therefore of the form 

H sdg = e d h d + e g n g - kQW ■ Q« - K (l - X 2 )Q^ ■ Q (4) , (11) 

where and are the quadrupole and hexadecapole operators, respectively, given 
by 

= a[s ] xd + d ] x s]W 

+ X [<$ x d]W - fa'[dt x~g + g^x d]^ + S&tf x ~g}^) , 

Qf = [Jx~g + tfxs]f- (12) 

The parameter x is chosen so that < x < 1 and a, a' = ±1 are two arbitrary 
phases leading to four possible realizations of SU(3) [5T]. For x — 1 the quadrupole 
operator is a generator of sdg-SXl(3) while for X — it becomes a generator of 
SO(15). Other parametrizations of the hexadecapole strength with the same effect 
could be taken in the hamiltonian (fill) , such as y/1 — x' 2 or (1 — \ X \); they differ 
little in actual results and the dependence (1 — x 2 ) is chosen here for its convenience. 
Note that with this parametrization the deformed limit has a vanishing hexadecapole 
interaction which is then only present in the 7-soft SO(15) limit. To do otherwise 
would require a hexadecapole operator taken from one of the other two strong coupling 
limits of sdg-IBM, SU(5) or SU(6). Since these have a questionable microscopic 
interpretation [52], this is not done here. 

With the help of the expression (JHJ) for the general hamiltonian of the sdg-IBM the 
classical limit of the special hamiltonian (ITT1) yields 

(H sdg ) _ e d /3* e g f3l (N - 1)« 



X 1 + + 1+% + Pl (l + /3 2 2 + /9 4 2 ) 2 

40* - ao'xfiyJtffiPi + ? vf [HPl + cos 3 7 2 

1 v 2( m$_Q4 I 5813/32 o2 , 3125 o4 , 18%/lQ M i a3a 

A 1^1372^2 "T 2058^2^4 ^~ 6174^4 ~T 343 [1±/J 2 /J4 

+ 10/%/3f]<7 / cos372 - eff^ [33/3|/3| + 35/3|] cos6 72 
4(iV-lKl- x 2 )/3| 



(13) 



(l+/? 2 2 +/?4 2 ) 2 ' 

where only the two-body parts of the ■ and ■ operators have been 
considered. We note first of all that a change of sign of the product <j<j' is equivalent 
with the change x ~^ ~X- We may thus take o~o~' = +1 and consider the parameter 
range — 1 < x < +1- Secondly, we remark that everywhere else a' occurs in combination 
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with cos 372 and, consequently, the change a 1 — > —a 1 is equivalent with the replacement 
72 — > rr/3 — 72. This means that without loss of generality the choice a = a' = +1 can 
be made as long as the range — 1 < x < +1 is considered. This can also be understood 
from the form of the QL operator in eq. (jT2l) from which it is clear that the change 
o~ — > — a is equivalent to x — ► — X- The two choices x = — 1 and x — +1 eacn correspond 
to a dynamical-symmetry limit of the SU(3) type which shall be referred to as SU_(3) 
and SU + (3), respectively. 

Choosing as the overall energy scale, we summarize the preceding discussion by 
stating that we seek to minimize (in (3%, /3 4 and 72) the following energy surface: 



V3(Aa,#4,72;»7,X>0 



ft r/3j 



x 4/3| - x^y/lftPt + f yfl [11/31 + 10/5 2 /3|] cos3 72 

1 _.2f 605 ,o4 , 5813/02/02 , 3125^4 , 18>/l0 r-i -i o3 a 
^ X \ 1372^2 "I" 2058^2/^4 6174*^4 ' 343 L 1J -^2M4 

+ 10/3 2 /3|] cos 3 72 - 6if4-[33/? 2 2 /? 4 2 + 35/31] cos6 72 ) 
4^7(1 - X 2 K 



(l + /?I + /? 4 2 ) 2 ' (14) 

where the notation rj = (N — l)«/e<j and r = e g /ed is introduced. The order parameters 
are fa, (3± and 72, and the control parameters are 77, \ and r. The subscript 3 in V3 
refers to the number of order parameters and is used to distinguish it from V*j introduced 
below. Physically meaningful values of the parameters imply 77 > 0. It is convenient 
to introduce, instead of 77 = £/4(l — £), the control parameter £ = 4r//(l + 4r/) such 
that £ G [0, 1] instead of 7/ G [0, +00 [. (The combination Arj instead of 77 is chosen in £ 
for later convenience.) Furthermore, since the g boson is higher in energy than the d 
boson we have r > 1. In the limit r — > 00 one recovers the sd-IBM. Again, it is more 
convenient to introduce, instead of r = 1/(1 — £), the control parameter £ = (r — l)/r 
such that £ G [0,1] instead of r G [l,oo[. The definition adopted here for 77 or £ is 
similar to the one of Iachello [6], so that we can easily compare with sd-IBM for ( = 1. 

The parameter space covered by the hamiltonian ( 1TT1) can be summarized with a 
generalization of the Casten triangle which becomes a prism, as shown in figure [TJ Also 
indicated on the figure are the different symmetry limits of the sd-IBM and sdg-IBM. 
The point (£, C) = (0, 1) is denoted as U(5), as should be, and the point (£, £) = (0,0), 
somewhat arbitrarily, as U(9), so that all vertices of the prism are labelled, facilitating 
the discussion of the structure of the phase diagram. Note that the U(5)®U(9) symmetry 
is valid for any combination of and e g (and hence for all ( G [0, 1]) as long as £ = 0. 
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Figure 1. The Casten prism associated with the hamiltonian (jlljl of the sdg-lBM. 
The different symmetry limits of the sd-IBM and sdg-lBM are indicated with black 
dots. The three control parameters £, x an d C as defined in the text, are indicated 
alongside the corresponding axes of the prism. 



5. Method 



The problem here is more complicated than in s<i-IBM because of an extra order 
parameter ((3^) and an extra control parameter (r). The necessary and sufficient 
conditions for the potential energy V to have a (Morse or normal) critical point in 

72*) are 



d(3 2 



0V 3 



d(3 4 



dV 3 



#72 



0. 



(15) 



(/3|./3I,7|) 



Furthermore, whether the critical point is a minimum, maximum or saddle point depends 
on whether the eigenvalues of the stability matrix in /3|, 72), 



H 3 ({3*, Pt,^ Z,X,0 



\ d*V 3 


d 2 V 3 


d 2 V 3 I 






<9/?2<97 2 


d 2 V 3 


d 2 V 3 


d 2 V 3 




dffi 


<9/?4<972 


d 2 V 3 


d 2 V 3 


d 2 V 3 


. <972<9/? 2 


d<y 2 df3i 


dil J 



(16) 



are all positive, all negative or both positive and negative. If, in addition to the 
condition (fl5l) . the determinant of the stability matrix vanishes, 

MPl PI 7 2 *; £, x, = det H 3 (p;,pt 72 *; £, X, = 0, (17) 

a so-called degenerate (or non-Morse) critical point is found [53]. This ensemble is of 
particular importance since it defines the points at which the energy surface changes in 
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structure and hence at which the system possibly undergoes a shape phase transition. 
The points that satisfy the condition (ITT)) define a surface in the control parameters 
space (£, x>0- Degenerate critical points shall be denoted as (£*,X*>C*)- O ur goal will 
be to determine the appearance of the degenerate critical surface in the Casten prism 
of figure [U 

The problem of determining the phase diagram for the Casten prism can be 
addressed by noting that solutions of the last of the equations in ( TT5|) satisfy one of 
the two equations 



sin 372 = 0, 



cos 372 



9/? 2 (9V2 X /?4 - 7V7)(ll/3f + 10/3|) 



(18) 



5^X^4(33/31 + 35/31) 

These two conditions define the two possible classes of solutions. The first class 
corresponds to critical points with axial symmetry, 72 = (prolate) or 72 = tt/3 (oblate). 
A further simplification is then obtained by noting that the off-diagonal elements of the 
stability matrix involving 72 vanish at 72 = or 72 = 7r/3, and hence the points in the 
degenerate critical set (£*,x*)C*) satisfy one of the two following equations: 



det 



d 2 V 3 

d 2 V 3 
d(3 4 d(3 2 



d 2 V, 



dp 2 dp 4 
'3 



d 2 V, 



0. 



d 2 V, 



0. 



(19) 



(/3 2 *A*>7 2 *) 



Furthermore, we note that cos 372 always occurs in Eq. (1141) in combination with an odd 
power (1 or 3) of p 2 . Since the replacement 72 — » tt/3 — 72 corresponds to a change in 
sign of cos 372 and no change in cos 672, both cases 72 = and 72 = tt/3 are covered 
by putting 72 = and allowing also negative values of p 2 . This further simplifies the 
potential energy surface (TT4")) which now only depends on the two order parameters (3 2 
and /?4, 



V 2 (P 2 ,Pi;ri,x,r) 



Pi 



+ 



r/3j 



l + Pl + Pi l + P 2 2 +P 2 4 (l + (3 2 + (3 2 y 
x 4/32 2 + x(f^ 2 3 + ? V / f^ + fv / f^ 



605 , 198V10 /J3/3 i 923,02/32 _i_ 180V10 
1372^2 ' 343 ^2M4 -T 343^2^4 343 



+ x 2 ( 

477(1 - x 2 )Pl 
(i + pl + piy 



P2PI + iif/^l) 



(20) 



In summary, critical points ((3%, 01) with axial symmetry are obtained from the 
simultaneous conditions 



dV2 
dp 2 



OP, 



0. 



(21) 
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where one assumes 72 = and allows negative values of (3%. The points (£*,x* 
the degenerate critical set should satisfy one of the two following equations: 



A 2 (/3 2 *,&*;£ )X ,C) = 0, 



d 2 V, 



Oil 



0. 



Ill 



(22) 



where A 2 (P2, PI; £, X-> is the determinant of the 2x2 stability matrix, 



W 2 (/3 2 *,/3 4 *;£,X,C) 



d 2 V 2 

d 2 Vo 



d 2 V 2 
dp 2 dp 4 
d 2 V 2 
dPl 



(23) 



dp 4 dp 2 

The second class corresponds to critical points without axial symmetry, obtained 
as solutions of the second equation in ( fl8l) . 

Before attacking the problem of the occurrence of shape phase transitions for the 
entire hamiltonian ffTTl) . it is instructive to study, as a function of r or (, the three legs 
of the triangle between the limits U(5) ®U(9), SU + (3) and SU_(3). This will determine 
the structure of the phase diagram on the faces of the prism of figure HJ In this analysis 
care should be taken to cover the two cases x > and x < 0, and study in particular 
the behaviour at x = 0, for instance in the SU_(3)-SO(15)-SU + (3) transition. 



6. Partial phase diagrams 

6.1. The U(5) ® U(9) to SO(15) transition 

The U(5) ® U(9)-SO(15) transition is obtained for x = and r = 1 (or ( = 0) while 77 
varies from to +00 (or £ from to 1). At the top side of the prism, for r = 00 (or 
C = 1), this reduces to the U(5)-SO(6) transition. We are interested here in the entire 
phase diagram in 77 and r (or £ and () which can be established by an expansion in p 2 
and Pi around (0,0). Since there is no dependence on 72 for x — 0, the analysis can be 
done starting from the potential ( 1201) . Up to fourth order we obtain 

V 2 (P 2 , A; V, X = 0, r) « (1 - 477)/?| + (r - Ar,)Pl - (1 - 8 V )p 4 2 

_ (1 _ i Qr] + r )02pz _ ( r _ 8^4, (24) 

involving even powers of p 2 and P^ only. The expansion (j24l) shows that the potential 



always exhibits a minimum, saddle point or maximum at P 2 = P% = 0. This is confirmed 
by computing the stability matrix at P 2 = PI = 0: 



^2(0,0;£,x = 0,C) 



2(1 -2Q 




2[1 





e(2 



01 



(25) 



(i-O(i-C) 

The zeros of the determinant of the stability matrix determine the degenerate critical 
points which separate the different possible shapes of the potential. They are £ = 1/2 
(or 77 = 1/4) and £ = 1/(2 — £) (or 77 = r/4). Consequently, we find that the potential 
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has (I) a minimum for £ < 1/2, (Ha) a saddle point for 1/2 < £ < 1/(2 — £) and (lib) 
a maximum for £ > 1/(2 — £), always at /3| = /3| = 0. 

In addition, for £ > 1/2, critical points in /3 2 occur for /3 2 7^ (i.e. quadrupole- 
deformed critical points). The dependence of /3 2 on £ and £ is obtained by solving 
dV 2 /d(3 2 = for & = 0, yielding 



/?* = ±V2£ - 1 = ± 1 



(26) 



4 + e d 

Hence the picture emerges that the spherical minimum /3| = /3| = at £ = 1/2 
continuously evolves into the two quadrupole-deformed minima (1261) . The phase 
transition is of second order. It should be emphasized that this result is independent of 
( (or r) and that the entire line £ = 1/2 exhibits a second-order phase transition when 
crossed. In particular, we recover the results of sd-IBM [6] which corresponds to ( = 1. 
One should verify the nature of the critical points by calculating the stability matrix at 
p; = ±V2£ - 1 and (5\ = 0: 

2£ - 1 



2(1 -£)£ 3 






C 



(27) 



/2£-l,0;£, X = 0,C) 

n 

Both matrix elements are positive definite and hence one always has a minimum. 

Finally, for £ > 1/(2 — (), additional critical points with /3| ^ occur 
hexadecapole-deformed critical points). The dependence of fi\ on £ and ( is obtained 
by solving dl^/dfi^ = for /3 2 = 0, yielding 



i.e. 



Pt = ±1 



'£(2-0-1 



i-£C 

So the spherical saddle point /3 2 = 



±4 



4 



(2? 



V 4 + e 9 

/?4 = at £ = 1/(2 — £) continuously evolves into the 
two deformed critical points ( 1281) . the nature of which follows from the stability matrix 
at f3* 2 = and # = ±^[£(2 - C) - W - £C) = 



^ 2 (0,± 1 



/£(2-0-l 



C(i 



i-£C 
-£C) 



ea-o 5 



o 



£,x = o,0 



o 



[£(2-C)-l](l-£0 J 



(29) 



2(1-£)£3(1-C) 4 

One matrix element is positive definite while the other is negative definite and hence 
both critical points are saddle points in this case. 

The phase diagram in the (£, () plane with \ = deduced in this way is shown 
in figure [2J The parameter space is divided in three areas with qualitatively different 
potentials in /3 2 an d /?4, characterized by the number and nature (minima, maxima 
or saddle) of Morse critical points. The potential surfaces in the regions indicated in 
figure [2] are characterized by 
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Figure 2. The phase diagram in the (£, £) plane of the Casten prism for x = 0. 
The lines or curves represent the locus of degenerate critical points that separate the 
different regions as discussed in the text. Symmetry limits of the sd-IBM and sdg-IBM 
are indicated with black dots. The entire £ axis corresponds to the U(5) <g> U(9) limit 
of sdg-IBM. 



• (I) a spherical minimum (/JJ = PI = 0); 

• (Ha) a spherical saddle point and two quadrupole-deformed minima 7^ 0, 

PI = 0); 

• (lib) a spherical maximum, two quadrupole-deformed minima and two 
hexadecapole-deformed saddle points = 0, PI ^ 0). 

The degenerate critical lines in the phase diagram correspond to bifurcation lines: at 
£ = 1/2 the spherical minimum splits into two quadrupole-deformed minima and at 
£ = 1/(2 — C) the spherical saddle point splits into two hexadecapole-deformed saddle 
points. 

6.2. The U(5) <g> [7/5/ to SU(3) transitions 

The U(5) ® U(9)-SU + (3) transition is obtained for x = +1 and r = 1 (or C = 0) while 
?7 varies from to +oo (or £ from to 1). A similar transition, from U(5) ® U(9) to 
SU_(3), is obtained for x — ~ 1- At the top side of the prism, for r = oo (or £ = 1), 
the U(5) ® U(9)-SU±(3) transitions almost reduce to U(5)-SU±(3), their analogues in 
sd-lBM. This correspondence is not exact since the coefficient of the [eft x d]ffi term in 
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Figure 3. The phase diagram in the (£,C) plane of the Casten prism for \ = —\. 
The lines or curves represent the locus of degenerate critical points with 72 = that 
separate the different regions as discussed in the text. Symmetry limits of the sd-IBM 
and sdg-IBM are indicated with black dots. The entire ( axis corresponds to the 
U(5) <g> U(9) limit of sdg-IBM. 

the quadrupole operator (Q3) is llvTo/28 « 1.24 while it is v^/2 » 1.32 in sd-IBM. 
From the expression of the potential surface (fT4l) we note that the change \ — > — x is 
equivalent to the simultaneous change ft — > —ft and 72 — > 7r/3 — 72. Hence all potential 
energy surfaces that can be realized for x < follow from those with \ > and only 
one case needs to be studied. We are interested here in the entire phase diagram in r\ 
and r (or £ and () and consider the case \ = — 1. 

An expansion up to third order around (fa. At) = (0,0) gives 

V(fo, ft, 72; ?7, X = -1, r) » (1 - 477)$ + r/3j - f^vP* cos 3 72 

- f yf^lft - f Vf^^ft 2 cos3 72 . (30) 

This shows that the spherical point Pa) = (0,0) is always a minimum in ft since 
r > and that it is a minimum in ft for £ < 1/2 (or 77 < 1/4). 

In general, critical points can only be obtained through a numerical solution of the 
equations (Tl5|) and (fT7|) . which, for a given x, can be achieved as follows [with reference 
to the classical limit (fl~3l expressed in units of n(N — 1)]. Since the hamiltonian is linear 
in and e g , it is straightforward to obtain as solutions of the first two equations in ffT5l) 
analytic expressions e<i(ft,ft) and e g (ft,ft). Furthermore, the solution of the third 
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equation in (1151) is also known in closed form, see (1181) . As a result, the "unknowns" 
€d, e g and 7 2 can be eliminated from f|T7|) to yield a single equation in /3 2 and (3^. This 
equation cannot be solved analytically but, for a given value of fa, and always for fixed 
given x, a numerical solution in /? 4 can be found, and from there 72, and e g are 
determined. Transforming from the (ed,e g , k(N — 1)) to the (£, () variables, one thus 
finds parametric curves in (£, Q space with /? 2 as parameter. Identical results are found 
by inverting the roles of f3 2 and /3 4 , in which case the parametric curves in (£, () space 
have /?4 as parameter. 

The solution for \ = — 1 already displays a surprising richness which is illustrated in 
figure [3] The different curves in the figure represent the locus of degenerate critical points 
with 72 = 0, that is, degenerate critical points that satisfy the first of the equations f|T8|) . 
The physical portion of the phase diagram with < £ < 1 and < ( < 1 displays 
several regions distinguished by qualitatively different potentials. We may characterize 
a potential by listing all its critical points together with the signs of the eigenvalues of 
the stability matrix which distinguishes minima from maxima from saddle points. Hence 
we introduce for each critical point a notation (Ai,A2,A3) with Aj = ±. A minimum, 
for example, is denoted as (+,+,+). A special situation arises for spherical critical 
points with (P^PD — (0,0) which can be characterized by the signs of the eigenvalues 
(Ai,A2) of the 2x2 stability matrix (I23I) . The different regions of figure [3] then have 
the following potentials: 

• (la) (+,+), (+,-,-), (-,-,-); 
. (Ib) (+,+), (+,+,-), (-,-,-); 

• (II) (+,+), (+,+,+), (+,+,-), (+,-,-), (-,-,-); 

• (Ilia) (+,-), (+,+,+), (+,+,-), (+,-,-), (-,-,-); 

. (Illb) (+,-), (+,+,+), (+,+,-), (+,-,-) 2 , (-,-, -) 2 . 

We emphasize that this is not the complete classification since it only concerns critical 
points with 72 = 0; an even more complex diagram is found if solutions from the second 
of the equations f|T8|) are included. In the study of phase transitions of quantum systems 
we are, however, primarily interested in the evolution of global and local minima of the 
potential. Therefore, we may for our purpose here ignore the evolution of saddle points 
and maxima, and concentrate on minima. From this point of view the phase diagram 
is divided into three regions characterized by potentials with (I) a spherical minimum, 
(II) a spherical and a deformed minimum, and (III) a deformed minimum. Also, if we 
confine our attention to the global and local minima of the potential only, no degenerate 
critical point with 72 7^ is found to occur in the physical region of the phase diagram 
with < £ < 1 and < £ < 1. From now on our analysis will be confined to the 
partition of the phase diagram into regions with potentials that differ in the properties 
of their minima. 

The phase diagram for the U(5) <S> U(9)-SO(15) transition of subsection 16.11 is 
characterized by two main regions (I) and (II) (see figure [2]) since (Ha) and (lib) are 
distinguished by differences in saddle points and maxima only. With respect to the result 
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Figure 4. The evolution of the deformation parameters at equilibrium, /3| and /3| , 
along three different paths in the plane with < £ < 1 and £ = 0, 0.5 and 

0.9, respectively. In each case the region of coexistence of a spherical and a deformed 
minimum is indicated by the grey area. 



of subsection [6J] we note that the U(5) ®U(9)-SU(3) transition has an extra region (II) 
which is characterized by the coexistence of a spherical and a deformed minimum. This 
result is consistent with what is found in the s<i-IBM and, moreover, we observe that, 
for realistic boson energies with e g > e^, the actual value of the g-boson energy has little 
influence on the size of the coexistence region. 

Although one does not have analytic expressions for the deformation parameters at 
equilibrium, /?| and (3%, these can be obtained numerically and typical results are shown 
in figure H] for three different paths in the (£, C) plane. The most striking evolution is 
that of the behaviour of (3± as a function of (. It is clear from the figure that a first- 
order transition occurs in (31 but the discontinuity in (3\ becomes smaller for £ — ► 1 (or 
e g -> oo). 

6.3. The SU^(3)-SO(15)-SU + (3) transition 

The SU_(3)-SO(15)-SU + (3) transition is obtained for 77 = +00 (or £ = 1) while \ 
varies from —1 to +1. In the limit rj — > +00 there is no dependence of the potential on 
either of the boson energies or e g , and hence the search for degenerate critical points 
is independent of (. It is easy to show that for x — the equations ffl"5l) are satisfied 
for P^ 2 + PI 2 = 1 an d arbitrary 7^. The potential is not only flat in 72 but in addition 
behaves like a Mexican hat in (/?2, At) around (0,0). Furthermore, for fi^ 2 + PI 2 = 1, 
the condition (fT7|) for a degenerate critical point is satisfied as well. Thus we find that 
all points with \ = and arbitrary ( are degenerate critical. On this line the potential 
is flat in 72 and in the SU_(3)-SO(15)-SU + (3) transition it tilts over from a prolate 
minimum with 7! = for % < to an oblate one with 72 = ir/3 for % > 0. No other 
degenerate critical points occur for values of x between —1 and +1. 
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7. The complete phase diagram 

We are now in a position to piece together the entire phase diagram for < £ < 1, 
— 1 < X < +1 an d < £ < 1. Since the change \ —>■ —\ is equivalent to the simultaneous 
change /?4 — > — /?4 and 72 — > tt/3 — 72, it is in fact sufficient to consider only positive 
or negative values of \'- the phase diagram is mirror-symmetric with respect to the 
U(5)-U(9)-SO(15)-SO(6) plane. In the following we consider the range — 1 < x < 0. 
An expansion of the potential up to third order, 

Ab 72; V, X, r) « (1 - 477)^1 + [r - 477(1 - X 2 )]/3f 

22 /5", y<r ,a3 „„„ 1 72 /2, /m o2o , 20 /5"„ 



+ f Vf^cos 37 2 + Uy/ixnffiPi + f V§X^Wcos3 72 , (31) 
or, alternatively, the expression for the stability matrix at /3| = /3| = 0, 

2(l-2£) 



^ 2 (0,0;e,x = 0,C) 



1-e 

Q 2{_l-£[2-C-X 2 (l-C)]} 



(32) 



(i-O(i-C) 

shows that at the spherical point the potential has (I) a minimum for £ < 1/2 (or 
7] < 1/4), (Ha) a saddle point for 1/2 < £ < l/[2-C-x 2 (l-C)] [or 1/4 < r] < r/(4-4x 2 )] 
and (lib) a maximum for £ > l/[2 — £ — x 2 (l — C)] [or 77 > r/(4 — 4% 2 )], always at 
P2 = PI = 0- The situation is similar to the transitions discussed in subsection I6.1[ 
with the only difference that now the separation line between the regions (Ha) and (lib) 
depends on the value of x while for the U(5) <S> U(9)-SO(15) transition we have x — 0. 
If we confine ourselves to an analysis of the properties of the minima of the potential, 
we may ignore the separation between (Ha) and (lib) and treat them as a single region. 

A locus of points of particular interest is defined by the condition that two 
eigenvalues of the stability matrix vanish. This is possible if there are at least two 
order parameters and at least three control parameters, as is the case presently, and 
leads to a catastrophe function of the type D±4 [53]. From the stability matrix (|32|) it 
is seen that this happens at /?| = (3% — for £ = 1/2 (or 77 = 1/4) and x 2 = C/(C — 1) 
(or x 2 — ( e d — e 9 )/ e d)- The locus of D ±i catastrophes is represented by a curve in 
the (£, x, C) space but it is entirely situated in the non-physical portion of this space 
(for e g < €d, below the Casten prism of figure [[]) with the exception of the single point 
£ = 1/2, x — C — which is half-way between U(9) and SO(15). 

We therefore come to the conclusion that nothing spectacular happens within the 
Casten prism which, when it comes to the properties of minima, is divided by two 
surfaces in to three regions (see figure [5]). The first surface is a plane obtained for 
£ = 1/2,— l<x<0 and < £ < 1. The second is a parametric surface with parameters 
ihiX) [° r (&> x)]j constructed with the technique explained in subsection 16.21 This 
situation closely parallels what is obtained in sd-lBM and, moreover, it is seen from 
figure [5] that the size of the coexistence region is not greatly influenced by the ratio of 
boson energies e g /ed (or (). We want to emphasize, however, that this is not a trivial 
conclusion but that this result depends on our choice of hamiltonian and of the domain 
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Figure 5. Phase diagram of the Casten prism associated with the hamiltonian (fTTj) 
of the sdg-lBM. The two planes represent the locus of degenerate critical points that 
separate the spherical and deformed regions, and enclose the region of coexistence. 
The different symmetry limits of the sd-IBM and sdg-lBM are indicated with black 
dots. 



available to the parameters £, x an d C which are guided by physics arguments. If we 
were to allow, for example, boson energies with e g < ej, a much more complex phase 
diagram would be uncovered which strongly departs from the one found in sd-IBM. 

8. Conclusion 

We summarize the three main assumptions at the basis of our geometric analysis of the 
sdg-lBM. 

(i) The general parametrization of a shape with quadrupole and hexadecapole 
deformation is reduced to one in terms of three variables (fa, fa, 72)- 

(ii) A simplified sdg hamiltonian is considered which is intermediate between four 
dynamical-symmetry limits of U(15), namely U(5) ® U(9), SU±(3) and SO(15). 

(iii) Parameters of this simplified sdg hamiltonian are restricted to physically acceptable 
values. 
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Under the above assumptions we find that the phase diagram is divided into three 
regions characterized by potentials with (I) a spherical minimum, (II) a spherical and a 
deformed minimum, and (III) a deformed minimum. The deformed minima are either 
prolate or oblate depending on the sign of the parameter \ m the quadrupole operator. 
No transition towards a stable triaxial shape is found. 
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